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It is shown in this paper that the charge-control concept can be conceived 
as a special form of the Linvill model for semiconductors. Instead of 
mathematical tools, charge-control models become equivalent circuits amen- 
able to ordinary network analysis techniques. In the simplest form, the 
charge-control equivalent circuit for the junction transistor is fully equiv- 
alent to the Linvill and the Beaufoy-Sparkes model. For all practical 
purposes, it is also equivalent to the Ebers-Moll model. 

The charge-control junction transistor equivalent circuit combines those 
features of the other models that are important for electrical engineering 
applications. It also permits the conversion between the three basic types 
of models. Because of its close relationship to the physical processes governing 
a device, it can readily be extended to higher-order phenomena. This is 
usually done by expressing a Linvill-type lumped model in terms of charge 
parameters. The charge-control equivalent circuit can be useful for modeling 
a variety of semiconductor devices. 

I. INTRODUCTION 

Three basic approaches are generally used to obtain descriptive large- 
signal models for transistors and diodes, the Ebers-Moll model, 1 the 
Linvill model 2 and the charge-control concept 3 after Beaufoy and 
Sparkes. 

The Ebers-Moll transistor model 1,4 is based on the idea of super- 
imposing a "normal" and an "inverse" transistor. Semiconductor 
junctions are represented by means of diodes and capacitors, whereas 
the properties of the transistor base are represented by frequency- 
dependent current sources. The Ebers-Moll transistor model is the 
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most popular of all transistor models since it lends itself most readily 
to simple "rule-of -thumb calculations." The current relations are 
described in the frequency domain, whereas the junction voltages are 
described as functions of current in the time domain, or, as in the 
original paper, only at dc. The model simulates only the effect which 
minority carrier storage exercises on the relations among the various 
device currents, but not the effect on current-voltage relations. Since 
the diode is a one-port device, no diode model of the Ebers-Moll type 
exists that could simulate carrier storage.* 

The Linvill model 2 ,5 ~ 13 is almost a direct representation of the con- 
tinuity and diffusion equations for semiconductor materials. It uses 
physical rather than circuit parameters and is superior to any other 
model when it comes to incorporating second-order physical effects 
or symbolizing new structures. 

The charge-control concept 3 ' 14 " 33 stands about halfway between 
physics and circuit considerations. It has proven in the past to be 
very useful for studying storage effects in diodes and transistors, but 
appeared to be entirely a mathematical tool. Certain equivalent circuits 
have been presented 14 ' 15 to illustrate charge control, but, as Linvill 
phrased it, "they have little more meaning than a symbolic model 
useful for the purposes of visualizing only." 

Hamilton, Lindholm and Narud compared the three models for the 
transistor in a well-written tutorial paper. 010 They discussed the 
approximations used in deriving each model from the same physical 
background. [See also Ref. 34] In contrast to this parallel treatment 
of the three models, the following study dwells on the interrelations 
and conversions between the various models. This is illustrated sym- 
bolically in Fig. 1. 

We may call the Linvill model a physical model, the Beaufoy-Sparkes 
charge-control model a mathematical model, and the Ebers-Moll model 
an electrical model. The link between the three models is accomplished 
through a modified approach to charge-control theory: instead of 
deriving, from device physics by means of integration, mathematical 
charge-control expressions, the charge-control concept can be treated 
entirely as an equivalent circuit tool. 27 The transistor model, for ex- 
ample, is in such a form readily comparable with, and convertible 
into the Linvill and the Ebers-Moll model, provided all of these models 
are at the same level of approximation. In its simplest form, the charge- 



* Diode models that simulate storage and use neither the charge control nor the 
Linvill concept are usually extensions of small-signal models towards incorporating 
certain nonlinear properties. 
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Fig. 1 —Principle of derivation of transistor models and their interrelations (heavy 
lines indicate main aspect of this paper; numbers refer to conversion equations in 
the text). 



control equivalent circuit model is fully equivalent with the standard 
form of the Beaufoy-Sparkcs charge-control model. But equivalency 
is usually lost, as extensions to higher-order approximations are made 
in each model. 

In this paper, we shall review the derivation of the above-mentioned 
types of models for diodes and transistors. This will be done with the 
help of a differential transmission line model. The equivalent circuit 
type charge-control concept will then be derived for diodes and tran- 
sistors. This will be followed by a discussion of higher-order approxi- 
mations, the inclusion of drift fields, and possible applications to 
other semiconductor devices. 



II. DIODE MODELS 

2.1 Mathematical Description 

As a starting point for our discussion it is assumed that the reader 
is familiar with the continuity and transport equations, describing 
current flow and carrier density in a semiconductor material. 
Continuity equations 
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— div j p (/) = e — — + e ~ U*J 

+div ;„(0 = e — ff + e U">J 

Transport equations 

j v (i) = evJ&pit) - eD p grad p(0 (2a) 

;.(0 = e/i»^i(0 + eZ) n grad n(/) . (2b) 

j p and ;'„ are the hole and electron current densities, respectively. 
p and n are the hole and electron carrier densities with p and n being 
their equilibrium values at a given temperature. E is the electric field 
intensity. D p and D„ are the hole and electron diffusion constants, 
and n„ and n„ are the respective carrier mobilities, e = +| e | is the 
value of the electronic charge. 

2.1.1 p-w Junction 

A p-n junction is described in a first-order approximation by the 
transport equation (2). The well- justified assumption is made that 
both ; P and j n are numerically small compared with the mutually 
opposing diffusion and drift currents. With the help of the Einstein 
relations 

D p = —» p (3a) 

D n = — M „ (3b) 

e 

and the appropriate boundary conditions one obtains the Boltzmann 
equations that express carrier densities as functions of the applied 
junction voltage v t 



^«t 



Pn(0, t) = p n0 exp [^ « ast (0 J (4a) 

n p (0, = n p0 exp I j^ t/„ t (0 J- ( 4b ) 

Here, p«(O,0 and n p (0,0 are the carrier densities on both sides of the 
junctions; p n0 and n p0 are the densities for v oxt = or, in other words, 
at points away from the junction, previously called p and n in (1). 
The definitions of these notations are illustrated in Fig. 2. 
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Fig. 2 — Carrier density distributions in the vicinity of a p-n junction and explan- 
ation of notations used. 

Iii terms of excess carrier densities, (4) transform into the following 
expressions 

P.,o...(0 = P»(0,t) - p nn = p n0 \exp<r^pv cxt (t)> - 1 (5a) 

n« B .„(0 = n»(0,0 - n p0 = n,o[exp|~« eJct (0| - 1 • (5b) 

Together with the reasonable approximation that the hole and electron 
currents pass through the junction unchanged,* (5) uniquely char- 
acterizes the junction. 

2.1.2 p and n Regions 

The following assumptions are implied in the analysis presented 
for a p-n diode : 

(i) The p-region is so heavily doped that the electron current can be 
neglected and appreciable carrier injection occurs only in the n-region. 
{ii) The problem is reduced to one-dimensional variations along the 
x axis. 

(iii) Drift fields are neglected. (Their inclusion will be briefly dis- 
cussed later in Section 7.3.4.) 

(iv) Space charge neutrality is assumed. 

* This is not quite true for silicon diodes at low forward currents and in the 
reverse direction where recombination in the space charge layer cannot be neglected. 
With respect to some of the diode properties, especially the current-versus-voltage 
relationship, the discrepancy can be accounted for by changing the exponent to 
ev„«/2kT.» 
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With these assumptions the continuity and transport equations reduce 
to 

d.v dt t v 

Jp(x,t) = -eD v — -- — V) 

We shall now express (6) and (7) in terms of the excess carrier den- 
sities p„„ „, which we shall denote for simplicity as p, i.e., 

P(x,t) = P„.„,..CM) = Pn(x,t) - P„n • 

Multiplying by the cross section A we obtain 

_ di P (x,0 _ eA &pM + Ci4 V&& (8) 

These are the two equations describing the n-region. 

2.2 Differential Diode equivalent circuits 

Equations (8) and (9) become transmission line equations if i v and 
p are taken as currents and voltages, respectively. (Mathematically, 
one may think of p as an analog voltage representing carrier density.) 
Fig. 3 illustrates the resulting r-g-c transmission line. 

The currents in the network branches are true currents but the 
voltages associated with the nodes are analog voltages. As a reminder, 
we have labeled the nodes with encircled "p's". The series and shunt 
elements are accordingly analog resistors, conductors and capacitors 
per unit length. 

If the diode is forward biased, the junction injects carriers into the 
n-region. They diffuse across the n-region gradually recombining until, 
a t x _> oo } a n hole current is converted into electron current. Fig. 4 
shows the carrier distribution across the n-region which is equal to 
the voltage distribution along the infinitely long r-g-c line. It can 
be derived easily from (8) and (9) that, under steady-state conditions, 
the shape of the charge distribution is proportional to 

exp(-.r/L„), 

where 

L D = \/D v t„ . (10) 
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Fig. 3 — Analog differential transmission line representation of diode model. (The 
bars indicate dimension per unit length.) 



L p is called the diffusion length. r v is the hole recombination time 
constant, or hole "lifetime". 

Since the analog voltage distribution on the capacitors of the r-g-c 
transmission line is identical with the physical charge density distribu- 
tion, and since many engineers have a much better feel for the charging 
and discharging processes of such a line than for the physical process, 
the r-g-c line representation may be quite helpful as an illustration 
of the carrier injection process. In early semiconductor work, such 
r-g-c transmission lines were frequently used. 20 ' 35 ' 36,37 No attempt was 
made, however, to attribute the physical meaning of carrier density 
to the network nodes; the junctions were represented by so-called 
/^-amplifiers. These amplifiers transform the internal voltage at x = 



P(x) 




Fig. 4 — Excess carrier distribution in diode n-region. 
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to the external voltage with the appropriate exponential relationship, 
while not transforming current at all. 

Linvill 2 ' 6 ' 8 ' 7 has introduced new symbols for the network elements 
which relate current to carrier density. These new notations avoid 
possible confusion between analog and physical circuit parameters, 
especially voltages, and hence enable us to combine current/carrier- 
density with current/voltage networks. Fig. 5 shows such a Linvill 
model in differential form. Again we have added bars over the letters 
as it was done with the f, g, and c in Fig. 3 to denote their dimensions 
as being "units per length". 

The symbols in the models are defined as follows : 



di v i(x,t) = —H c dxp(x,t) 
di v2 (x,t) = -S dx ^^ 



where 



a/H 



dp(x,£) = -(1/H d )dxi p (x,t), 

H c = combinance per length = eA/r p 
S = storance per length = eA 

1 



diffusance 



per length = ]/eAD v 



(lla) 
(lib) 
(lie) 

(12a) 
(12b) 

(12c) 




= ?4 + ^ 



h— • 



3 



Fig. 5 — Differential Linvill diode model. Note that, in consistency with common 
transmission line notations, the reciprocal of diffusance must be used. 
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This model can be extended to include majority carriers, drift fields 
etc. The reader is referred to the literature. 2,5-13 

2.3 Integrated Diode Models 

2.3.1 Mathematical Integration 

In order to arrive at an expression for the external diode current 
from the continuity equation (8) we can integrate this expression with 
respect to x. Choosing x = and x = oo as the limits of integration, 
we obtain 

r di p {xj) , r . d P (x,t) , i r , , N , 

-;„ V dx = I eA m dx + 7 9 1 eApM dx - w 

The third integral represents the total charge in the bulk material. 
With the appropriate boundary conditions i p (0) = i, i n (0) = 0, 
4(°°) = 0, *n(°°) = i, the well-known charge-control equation 3 can 
readily be obtained as 



at r„ 



(14) 



To obtain (14) from (13) the assumption must be made that A and 
r p are constant. Note that no approximations or restrictions to specific 
charge distributions are implied in (14). (They must be made, however, 
when relating the current to the junction voltage.) 

2.3.2 Lumped Linvill Diode Model 

The crudest approximation to the distributed Linvill model of Fig. 5 
is to replace the "line" by just one storance and one combinance 2,7 
as shown in Fig. 6. These two elements are obtained by summing, i.e., 




v(U 



v(t)=^ T ln 



Fig. 6 — Lumped Linvill diode model showing single-pole approximation for 
minority carrier storage. Chosen values: Ax = L p , p(t) = p(0,t). 
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integrating, all differential storances and combinances from x = 
to some value Ax. The value of Ax is usually chosen to equal the diffusion 
length L p . This may seem arbitrary, 12 but has no effect on the terminal 
properties of the first-order model, as long as p(t) is chosen such as 
to maintain the same amount of total charge. 

The values of the circuit elements follow from (11), (12), and (5a) as 



H, = H r Ax = H C L B = 



eAL„ 



S = S Ax = SL n = eAL B 



■w-W + 'S? 



where A is an abbreviated notation, used hereafter for 



X " kT 



(15a) 
(15b) 

(15c) 

(16) 



The meaning of such lumping with respect to the carrier distribution 
is illustrated in Figs. 7 and 8. The solid lines in Fig. 7 present the 
actual carrier distribution in a switching example in which a current 
pulse is assumed. As required by the transport equation, the slope 
at x = is, at any time, proportional to the current. Under steady- 
state conditions, an exponential distribution is obtained. To assume 
such exponential distributions at any instant of time (dashed lines 
in Fig. 7) represents a simplifying assumption. The corresponding 



p(x,t) 



p(x,t) 





(a) 

Fig. 7 — Illustration of the (a) charging and (b) discharging process in the neutral 
bulk material. The applied signals are assumed to be forward and reverse current 

pulses. The solid lines represent the actual shape for current pulse drive; the dashed 
lines represent exponential model approximations. 
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TOTAL CHARGE =q. (t) = eAI_pp (0,t) 

(a) (b) 

Fig. 8 — Exponential and corresponding lumped distribution of excess minority 
carriers in the bulk material of a diode, (a) Illustration of the choice of lumping 
length, (b) Time variation for m = 1 ; m = 1 is generally preferred in the Linvill 
model, and is irrelevant in charge models or circuit applications of the Linvill model. 



errors are negligible in all those applications where the switching times 
are large compared with the carrier redistribution times (= diffusion 
times r a and r b in Fig. 12). 

In the lumped Linvill model, it is assumed that the carrier density is, 
at any instant of time, constant from x = to x = L„ , and that it is 
for all x > L v . Any information on the distribution of the charge, 
especially of the slope at x ™ 0, as expressed in the transport equation, 
has been lost since all series elements (diffusances) are neglected. The 
only parameter of importance left is the total number of minority 
carriers and hence, the total charge. The approximation used is therefore 
equivalent to the dashed line exponential distribution in Fig. 7. 

As mentioned above and illustrated in Fig. 8(a), the length Ax over 
which p is nonzero, is most conveniently chosen to equal L p . But it 
is permissible to choose Ax ^ L p if the constant value p(x) is recognized 
to be different from p(0,t) ; for Ax = mL p , we must choose p(t) = p(0,t)/m 
such as to yield the same total charge 

q = eAL p p(0,t). (17) 

Fig. 8(b) shows, for m = 1, the time variation of the carrier distribu- 
tions for the lumped model (solid lines) and the exponential distribution 
(dashed lines). 

As the external voltage v(t) varies, the carrier density p(0,t) changes 
accordingly. The relation between v(t) and p(0,t) has been given above 
in (loc). We shall see below that the approximation made in the lumping 
process, as discussed above, effects only v(t) but not the current. Little 
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or no error is made whenever, and as long as the external driving 
source impedance is large. 

Since, in this and any other lumped Linvill model, all circuit param- 
eters are functions of total charges in the various sections of the semi- 
conductor we can transform the model of Fig. 6 into a charge-controlled 
model, in which all carrier densities are replaced by charges. This will 
be shown in Section 2.4.1. 

2.4 Charge-Control Diode Model 

Without invoking any of the approximations introduced in Section 
2.3.2 and illustrated in Figs. 7 and 8, the diode is completely described 
by (14) and (5a). These two equations are the basis for the classic 
charge-control theory after Beaufoy-Sparkes as applied to diodes. 
Throughout the charge-control literature, only the current appears as 
a function of the total charge but not the voltage. If we want to relate 
the junction voltage to the charge rather than to p n (0,t) as it was 
done in (5a), we must make some approximation: The simplest possible 
approximation is the assumption that p(Q,t) is proportional to q(f). 
This is, for example, satisfied if the shape f(x) of the carrier distribution 
never changes, i.e., if the carriers redistribute themselves instanta- 
neously. p(x,t) is then of the form 

p(x,t) = j(x)g(t). 

The shape of /(re) does not matter as long as the integral / f(x) dx 
yields the proper proportionality constant. Examples of this are the 
exponential distribution or the lumped distribution (with any arbitrary 
value of m) in Fig. 8(a). This shows the equivalency between the 
postulate of instantaneous carrier redistribution in classic charge- 
control theory and carrier density lumping in the first order Linvill 

model. 

If we now want to establish an equivalent charge-control circuit 
we must first represent (14) by corresponding circuit symbols. This is 
done in the n-region part of Fig. 9. S is the store originally introduced 
by Beaufoy and Sparkes. 11,12 To account for directionality, we have 
added a vertical bar to the store symbol in the manner of the standard 
diode symbol. The properties of S, as defined in this paper, are: 

(i) charge stored = q(t) 
(ii) current in direction indicated by arrow -^ ( ] 

; Ml 

1 " dt 
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JUNCTION n-REGION 



V(t) 
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PROPERTIES OF S: iT) L(t) = 



dq(t) 
dt 



OR I(s) = sQ(s)-Q(o + ) 

(Z) VOLTAGE = O 

Fig. 9 — Complete first-order charge-control equivalent circuit (this circuit is a 
charge representation of Fig. 6). 



(m) voltage across store = 0. 

S is often interpreted as an infinite capacitor for which i = dq/dt = 
d(Cv)/dt = finite, but C — > °o and v — > 0. 

It follows from (15c) and (17) that the junction voltage is of the form 

v(t) = £ In [1 + Kq(t)], 

where K is a proportionality factor. If we denote the steady-state 
reverse current (flowing through the diode when v(t) is very large 
and negative) by I s we can evaluate the constant: For v — > — qo we 
obtain 

KQ = -1 
and from Fig. 9 

Hence, 
and thus, 



-Is = Qh 



K = 



S T r 



■« ■ I "■ b + ft]- 



(18) 
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Under steady-state conditions where Q/t p = I, this equation becomes 
the well-known diode equation. 

For the possibilities of incorporating the junction capacitances see 
the discussion in Section 7.3.2. 

2.4.1 Derivation of the Charge-Control Model from the Lumped Linvill 
Model 

It represents not merely an additional proof of equivalency but also 
a good preparation for the derivation of more complex models, if we 
show 27 that we can derive the charge-control model from the Linvill 
model. A somewhat related modification of the Linvill model was more 
recently proposed by Beddoes. 12 To this end we calculate the currents 
through the elements H c and S in Fig. 6 : 

i„£t) = H r p(0J) (19a) 

is(0 = S^- (19b) 

Substitution of the values for H c , S, and p(0,t) from (15a), (15b), 
and (17) yields 

W« = ^pM = f (20a) 

isi0 = eAL /MA = «f. (20b) 

This result shows that the current source in the charge-control 
model of Fig. 9 represents the current i Be through the combinance, 
and that the store S represents the current i s through the storance. 

To find the expression for the junction, we can express p(0,t) in 
terms of q{t) by means of (17). p„ can again be obtained from the case, 
where V — > — », and where p(0,t) = P(0) = — p„ : 



/*Jr— - -Is = eA ^P(0) 



eAL„ 

~ P„n 






Thus, we find 

p(Q, t) = _q(Jl J Jsj, = q(t) 
p, lo eALj,/ eAL p I s t p 

With this we can make the transition from (15c) to (18). 



(21) 



2.4.2 Evaluation o/ the Charge-Control Model 

The charge-control model is completely equivalent with the lumped 
Linvill model in Fig. 6; in fact, it may be considered a circuit oriented 
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form of the Linvill model. In almost all instances 7 '' l where the Linvill 
model is being used for circuit applications the conversion of carrier 
density into charge must be made anyhow. The charge-control equiva- 
lent circuit in Fig. 9 uses current and voltage sources plus one lesser 
known circuit element described by the simple relations 

v(l) = (22a) 

m = d -f (22b) 

or, in Laplace notation 

I(s) = sQ(s) - Q(0 + ). (22c) 

Ordinary circuit analysis techniques can be used in working with 
the model. No restriction exists with respect to the external waveforms. 
Q appears as an additional circuit parameter with additional com- 
plexity comparable to that of an additional branch current. From a 
topological viewpoint it is a branch current. This is the price to be 
paid for inclusion of the first-order dynamic storage properties. 

Junction and n-region are clearly separated in the model. Thus, little 
difficulty should arise in adding junction capacitors (dashed in Fig. 9), 
series path resistors, and leakage resistors, provided, physical knowledge 
of such effects exist. 

2.4.3 Charge-Control Model for Short-Base Diodes 

Diodes with extremely short bases do not show the exponential 
minority carrier distribution represented in Fig. 4, but rather a prac- 
tically linear fall-off (like in a transistor base except that the collector 
is now a nonrectifying contact). With reference to Figs. 3 or 5, this 
means that the distributed "transmission line" is so short that the 
effect of the series diffusances H d dominates over that of the shunt 
combinances H c . The metallic contact behaves like a short circuit 
at the end of the line. 

The analogy with the r-g-c line of Fig. 3 may help the reader visualize 
the difference between the long base and the short base diode: The 
first-order approximation for the infinitely long line with respect to 
currents and input voltage is the parallel connection of the shunt 
resistor and the shunt capacitor; the first-order approximation for a 
very short line is the parallel connection of the series resistor and the 
shunt capacitor. In terms of the Linvill model, the short base diode 
model is obtained by replacing H c in Fig. 6 by H d = eADJL and 
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L v by L. (Note that H d increases as L becomes small.) In the charge 
control model, of Fig. 9, the term t v which represents the recombination 
time constant for the long base diode, now becomes the diffusion time 
constant. The new value t' v can be derived most easily from the Linvill 
model as follows: 



Q 



Q 



**, '" p(0)H d /W p(0)eAD p /W 2D, 



(23) 



Apart from this numerical change, the model in Fig. 9 for the normal 
diode is equally valid for the short base diode. 

2.4.4 Piecewise Linear Charge-Control Diode Model 

For many practical purposes the logarithmic voltage source relation 
can be approximated by a switch as illustrated in Fig. 10. The switch 
opens when q becomes negative and closes when q is able to charge 
up to q > 0. A threshold voltage V ih is connected in series with the 
forward path. If desired, the slope of the logarithmic curve 



dV dV 



dl d(q/r) \q Xg„ 



1 

X/ a 



x/, 



(24) 



v=f(q) 




Vth- 



v=f (q) 



-^"hSr) 



"SLOPE = -r-= 

Xq. 



q/r 




(a) 



(b) 



Fig. 10 — Piecewise linear approximation for the semiconductor junction, (a) 
Theoretical logarithmic curve, (b) Approximated curve (the dashed lines indicate the 
completion of the diode model). 
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can be added as a resistor, where I av is an average current, which may, 
in long hand calculations, be assumed to be %I mn * . The saturation 
current I s must now be represented by an external current source. 

2.4.5 Application of the Model 

The above discussion of diode models serves two purposes: First, 
they form a basic understanding for deriving transistor models. Secondly, 
the diode models can be very useful in simulating dynamic effects due 
to carrier storage in diodes. 

With the piecewise linear junction approximation of Fig. 10 applied 
to the charge-control model in Fig. 9, storage time equations can be 
derived easily using Laplace transform concepts. The model has proven 
to be very useful in the analysis of step-recovery diode circuits. In the 
piecewise linear form, it can be handled without a computer, whereas, 
for the more complex models with various parasitics added, computers 
soon become mandatory. 

Switching times for step-recovery diodes are derived in Appendix A.l 
as an example of the use of the charge-control model. The equations 
obtained have been found by many authors to agree well with actual 
measurements. The normalized storage time for recovery from an 
infinite ON-pulse according to (97) is plotted in curve a of Fig. 11 as 
a function of the reverse-to-forward current ratio according to the 
relation 

T = rln(l +^)- (25) 

When applying this result to an ordinary diode with homogeneous 
doping profile, one must be aware of the implied approximations: 
(i) The single-section approximation in the model does not affect any 
mutual relationships between currents and charges, but represents 
approximations with respect to the junction voltage. As the amount 
of stored charge is reduced considerably in the diode, the junction 
voltage decreases noticeably, (it) As the carrier density near the junction 
becomes extremely small, the voltage reverses sign and the diode 
impedance, at some point, becomes comparable with the external 
source impedance. The ideal current source assumed in (97) ceases 
to exist, and instead of the step-recovery, as given by the model, a 
long tail in the current response results. 

From either one of the two differential models in Figs. 3 or 5, we can 
calculate the time in which the carrier density at x = 0, and hence 
the junction voltage, reaches zero. Such a calculation yields the relation 
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Fig. 11— Comparison of diode storage times as functions of the driving ratios, 
(a) Single lump model; T = time when charge is fully depleted, (b) Differential 
model ; T = time when excess carrier density at x = reaches zero. 



originally derived by Lax and Neustadter 1 



(26) 



This relationship is illustrated in Fig. 11(b). Since curve (b) represents 
only the storage phase but not the very long tail of the recovery, the 
values are much smaller than those in curve (a) in which some sort 
of "effective total recovery" is represented. The difference is most 
remarkable at strong relative reverse drives where the carrier distribu- 
tions on the lines differ most from the steady-state distributions. 

If it becomes necessary to incorporate the tail of the recovery into 
a lumped diode model, the double 7r-extension described below may 
prove adequate for most applications. 
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2.4.6 Higher-Order Approximations 

Bearing in mind how the model originated as an approximation to 
the differential transmission line or as just another form of the lumped 
Linvill model we can now understand how higher-order approximations 
are to be obtained. 

Fig. 12 shows the example of a ^-approximation for a diode. The 




CLp = ^Lp 



(a) 



v(t) 



id 



aWn d,(t) *£•>([) q 2 (t) 

r - T r 




WHERE L d (t) = 



q.(t) 



q 2 (t) q,(t) q. 2 (t) 



acLpVDp bcL p 2/D p 



-)=>(.^') 



T v (FROM DC CONSIDERATIONS) = 



'■(■k-0+- 



CHARGE CONSERVATION CONSTRAINT ON a,b,C: 

cbO-a)L p 2 



(a + b-i)D p T 2 

WHERE C IS MOST APPROPRIATELY CHOSEN TO BE C = 

(b) 



a+b 

2 



Fig. 12 — Higher-order, Tr-Approximation of diode charge-control model, (a) 
Charge approximation, (b) Corresponding model. 
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charge is broken up into two parts q x and g 2 . The diffusance between 
the two stores controls the redistribution of the charge. Such a structure 
provides a better representation of the junction at the higher frequencies 
or at higher speeds than the model of Fig. 9, since the junction voltage 
is now a function of only that part of the total charge which is close 
to the junction. The model simulates recovery tails. It also permits 
the simulation of variations in recombination time along the z-axis. 
Fig. 12 assumes two different recombination times t x and r 2 . The 
i = /(g) relation then becomes 

. = d(g, + gj + fii + 22 (27) 

dt T\ r z 

[which reduces to (14), if one assumes r x = rj. 

Three additional time constants r„, t 6 and t, appear in Fig. 12. 
They depend on the choice of the sections aL p and bL p over which the 
shunt elements are integrated and on the choice of the section cL p over 
which the diffusances are integrated. The three degrees of freedom 
reduce to one, however, if one considers that (i) the total charge must 
be conserved by the lumped approximation, and (it) in a multisec- 
tional approximation the diffusances are most appropriately lumped 
over sections cL v which extend between the centers of the charge 
sections. The corresponding relations are given in Fig. 12; derivations 
have been omitted. 

III. LARGE-SIGNAL TRANSISTOR MODELS 

In complete analogy to the diode models, we shall now compare 
the various junction transistor models and establish the charge-control 
model in the form of equivalent circuits. The Ebers-Moll concept, 
which was found not to be applicable to dynamic diode description, 
will now enter the "competition". 

In order to dwell on the philosophies underlying each concept we 
shall, at first, limit ourselves to diffusion type junction transistors, 
neglecting again drift currents and secondary effects such as base-width 
modulation. All derivations will be carried out for pnp transistors; 
but, of course, everything will be correspondingly valid for npn tran- 
sistors. 

3.1 Differential Transistor Model 

The most rigorous of all the equivalent circuits describing a junction 
transistor, as defined by (5), (8), and (9), is the differential model shown 
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E r e ' 




1 = 
v e (t)=f[p(o.tj] 



i=W 

- v c (t) = f [p(W,tj] 



Fig. 13 — Differential Linvill model for the transistor with drift fields neglected, 
pnp version shown. 

in Fig. 13. Linvill notations comparable to the diode model in Fig. 5 
were chosen. (If so wanted, the model could also be drawn with the 
notations used in Fig. 3 resulting in an r-g-c line and two .^-amplifiers 
at both ends.) 

The base section of the transistor model is only a very short "trans- 
mission" line when compared with the "infinitely long" diode n-region 
of the normal diode. Instead of 100 percent recombination, as found 
in the diode, the transistor must have as little recombination as possible 
in order to achieve high gain. Fig. 14(a) shows a steady-state charge 
distribution under normal forward operation, and Fig. 14(b) shows 
the distribution for the case where both junctions are emitting, i.e., 




(a) 




Fig. 14 — Excess minority carrier distribution in the transistor base under (a) 
normal and (b) saturated operation. 
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in saturation. Under normal operation the collector acts as a "charge 
short circuit" for the line. For high-gain units, the slope is almost 
a straight line; at x = it is proportional to I E and at x = W propor- 
tional to I c • 

The general case is that of Fig. 14(b) where both junctions are 
emitting and p(W) ^ 0. Any section of the base region can be de- 
scribed analogously to a four-pole using the definitions given in Fig. 15. 
Note that there are no nonlinearities in the base section. 

I B (s) = A„P,(«) + A 12 P 2 (s) (28a) 

Io(s) = A 21 P l (s) + A 22 P-M- (28b) 

By using complete analogy to standard transmission line theory, it 
can be shown that with the use of (10) and (11) one obtains for a 
homogeneous section Ax 



J E (s) = — ^ coth 7 Ax — — jj-* cosech 7 Ax 



lc(s) = 



z 



cosech 7 Ax — 



Z 

PM 



coth 7 Ax, 



(29a) 
(29b) 



where 



Z = -r 

eA 



1 -{-st 



(30) 



(31) 



In the general case, the base is not homogeneous, which means that 
Z and 7 will vary along the line. 

The junctions are described by the time relations 

Pi(0,0 = P»o[exp {\vXt)\ - 1] (32a) 

p 2 0M = P«o[exp \\v c (t)\ - I]. (32b) 




Fig. 15 — Symbols and polarity conventions denning the four-pole description of 
the transistor base. 
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v e (t)=-^ln(i + 



v c (U=-^ln(.+ 



Pno 
P 2 W 



_ L|(t) eA,w, 
H c- - "pot) - 



"*" Pa(t) 



i 5 (t) _ eA 2 w 2 
r 2 



L 4 (t) 

S2= dMt)7dt =eA2W2 



H d = 



L 3 (t) eADr 



P,lt)-p 2 (t) w 
Fig. 16 — First-order lumped Linvill transistor model. 



IV. THE LUMPED LINVILL 7T-MODEL 

In the literature, any lumped approximation to the transmission 
line model presented in Fig. 13 is referred to as a "lumped Linvill 
model" or simply "lumped model". The most common form is the 
7r-model. With respect to its current properties and steady-state voltages, 
this form will prove to be equivalent to the commonly known "Ebers- 
Moll model", if both models are taken to be in the form of first-order 
approximations. 

By integrating all differential diffusances over a length Ax = W, 
all emitter sided differential combi nances and storances over a length 
Ax = Wi , and all collector sided differential combinances and storances 
over a length Ax = W 2 (= W — IFO one obtains the circuit shown 
in Fig. 16. Nonsymmetry has been taken into account by using different 
recombination times r, and t 2 and different cross-sectional areas Aj 
and A 2 on the two sides. Note that the latter represents an extension 
from the one-dimensional carrier flow and as such an example for the 
reduction of multidimensional effects to a one-dimensional model. 
Area A is some average cross section effective for the diffusion process. 
H d is the diffusance, the H c 's are the combinances, and the S's are 
the two storances. 
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The four-pole equations describing the transistor base are obtained as 



H£ft« 



r /n eAD v n ,J, . AJVWi,. . ' 

7 *(«) - if^^L 1 + ~4^at (1 + st,) . 

- tf d P,(s)[l + §^ (l + s |j-)] - W,(«) (33a) 

- JTA(t) - ffA(.)[l + §* (l + s J^) ] • (33b) 

The junctions are described as 

Pi(0 = P(0,0 = ?Uexp {Xw.(01 - 1] ( 34a ) 

p 2 (0 = p(PT,Q = p„ [exp \\v c (t)} ~ 1]. (34b) 

A constraint has to be satisfied: Under equilibrium conditions, the 
total charge in the base must equal that in the two sections, i.e., 

eA l W l P l + eA 2 W 2 P 2 = eA f P(x) dx W leAW[P(0) + P(W)]. (35) 

Jo 

The approximation holds for high-gain units. For this case the base 
volume sections are equal, i.e., AJVi = A 2 W 2 = \AW. For low-gain 
units (34) must be modified: The terms p v {t) or p 2 (t), or both, must 
be replaced by p,(i)/m, and p 2 (t)/m 2 , respectively, whereby the m's 
are constants > 1, similar to m in Fig. 8. 

Equation (33) represents one of several possible approximations to 
(29) with the additional property of nonsymmetry being added. 

Higher-order approximations of a lumped linear model are obtained 
by representing the base of width W by more than the two sections 
W l and W 2 . 

V. THE EBERS-MOLL TRANSISTOR MODEL 

The focal point of the Ebers-Moll model is the two-port description 
of the base. Such a description has been given in (28) and (29), and 
is permissible because of the linearity which exists between currents 
and carrier densities in the base. Nonlinearity exists, however, in the 
relationship between carrier densities and external voltages according 
to (32). Since linearity allows the use of the superposition principle, 
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the total current can be conceived as consisting of the superimposed 
contributions of the currents injected by the two junctions. 

When put into the form of an equivalent circuit, the Ebers-Moll 
model shows the superposition of a normal transistor (subscript N) 
and an inverse transistor (subscript /). In Fig. 17(a) the lower diode 
and current source represent the normal transistor and the upper 
elements represent the inverse transistor. Each junction is represented 
by a diode, a fraction of the diode current is collected by the other 
electrode. The ratios of collected currents to emitted currents are called 
a N and a r for normal and inverse operation, respectively. The general 
frequency behavior of the a's can be calculated for a homogeneous 
base from (29), (30), and (31) as 



OjIcpCsD 



/ev ( 
LcfK.)=I C fo| exp (- 



Eo VW 



lEF(t) = I EF0 [EXp(^Uj 



'e(t) \1 Sr b ' «nIef(s) 



BN 




' VW — °c 



(a) 



Eo WV 



-a r I c (s) L c '(t)=I C o[EXp(^|^-i)] 



ne 



LE 



W^-HT-)] 



■M * 



i-c' 



Ctc 



ne- 






ev e (t) \"| >r b ' «nIeCs) 



. , a NO , . Q IO 

WHERE Ot N (S)= o— Ot(S)- «. 

l + rr — 1+7- 2 — 

"on w qI 



(b) 



Fig. 17 — The two forms of the Ebers-Moll transistor model: (a) direct representa- 
tion of the idea of superimposing a normal and an inverse transistor, (b) collecting 
current sources as functions of the electrode currents. The junction saturation currents 
in (b) are identical with the open-electrode diode saturation currents. 
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cosech 7 Ax 1 



<8) = 



W«) 1 

^input(s) Jp„ n i p ut=0 

9.n t 

(36) 



coth 7 Ax cosh 7 Ax 
I 2Z) D r 



7 2 A.-r 2 2D p t + A.r 2 + sr Ax 2 
+ 2 

But symmetry does not exist in a practical transistor. The constants 
in (36) are therefore, different for a N and for a T . Equation (36) can 
be rewritten under this consideration in the well-known form 

•« - rrfe (37a) 

-<«> - r+fe;- (37b) 

The relations between the constants in (37) and the physical param- 
eters (corresponding to the constants in (36) modified for the non- 
symmetrical case) will be derived in Section 5.1. 

On account of their nonlinearity, the junction diodes must be de- 
scribed in the time domain. In their original paper, Ebers and Moll 
defined only a dc relationship between voltages and currents. This 
would restrict the use of their model to piecewise linear analysis. 
But the model can be made more general 118 by postulating that the 
v = j{i) relation be valid at all times, as indicated in Fig. 17. 

In either case, an important property of the semiconductor junction 
is lost: Voltages and currents appear as being directly related instead 
of being related indirectly through current density or charge. This can 
best be illustrated by an example. If a forward current through a junction 
is suddenly replaced by a reverse current the voltage actually does not 
reverse sign until the excess carrier density at the junction is reduced 
to zero. According to the Ebers-Moll model, voltage and current always 
change polarity together. As mentioned before, it is for this reason that 
for a diode, no dynamic model of the Ebers-Moll type exists that 
would represent charge storage effects. In addition to this shortcoming, 
the feature of mixed time and frequency domain characterization is 
undersirable if the model is to be used in its nonlinear form, say on 
a computer. 

The Ebers-Moll model was originally presented in a form, shown 
in Fig. 17(b), which differs slightly from that in Fig. 17(a). Both ver- 
sions have been used throughout the literature over the past years 
and very few authors 39 40 have clearly pointed out the difference between 



CHAROE-CONTROL CONCEPT 549 

them. In Fig. 17(b) the collecting currents are a times as large as the 
total emitter and collector currents, respectively. A simple calculation 
shows that the two versions are formally equivalent, if the relations 

'"« - 1 - Sbo (38a) 

and 

^"^ = 1 ^i M ( 38b ) 

1 — a N (s)ai(s) 

are satisfied. A glance at the equations for the voltage sources in Fig. 17 
reveals that the two versions could not be completely equivalent, 
unless either I EP0 and I CF0 or I E a and I C o would be considered frequency 
dependent. Due to the approximative nature of both models, this is 
normally not done. 

From both Fig. 17(a) and (b) the respective four-pole equations, 
on which the model is based, can readily be derived in terms of elec- 
trical parameters: 

I*® = I EF (s) ~ a,(s)I CF (s) = r f ~ a !^j s) (39a) 

I c (s) = a N (s)IMs) ~ Icris) = a i S) ^ ( ^f - (39b) 

After substituting the expressions for the junctions one obtains for 
the steady-state case the well-known Ebers-Moll equations 

Ie = r - [ ^ [exp (XV.) - 1] - aiJco [exp (\F f ) - 1] (40a) 

Ic = ; amlE0 [exp (\V.) - 1] - r -^ [exp (\V r ) - 1]. (40b) 

5.1 Convparison Between the Ebers-Moll and the Linvill Model 

Comparing (40) with (33) and (34) for the steady-state solution 
leads to the following relations: 

eADjp n0 _ asolEo _ arolco ,*-.} 

IF 1 — a m a I0 1 — aAT «/o 

A corresponding comparison for the ac case would yield the same 
expression as in (41), except that a, vo and a I0 would have to be replaced 
by their frequency dependent forms. Since the left side term of (41) 
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is frequency independent, no rigorous equality exists between the 
Linvill model and any of the two versions of the Ebers-Moll model 
under ac conditions. In the Ebers-Moll model, the junction voltage 
is a function of the total diode current [being different in the two 
versions of Fig. 17(a) and (b)]; in the Linvill model it is only a function 
of the resistive component of the diode current in Fig. 17(a); this 
component equals the current through the combinance which is propor- 
tional to the carrier density p. It can be shown that the correct solution 
in which the junction voltage is a function of the carrier density directly 
at the junction, lies between these two cases but much closer to the 
lumped Linvill simulation. The discrepancy, mentioned here, affects 
only the junction voltages and does not appear in many analyses 
that use piecewise linearity. 

a N (s) and a r (s) can be expressed in terms of the physical parameters 
by comparing (39) and (33) separately for the normal operation (I CF = 0) 
and for the inverse operation (I EF = 0). Subsequent conversion of 
the a's into /S's yields 

Mi > = r^Ai) = t^t - i + st, (42) 

a t v «M - P" - At 2 DJA 2 WW 2 , 

M) " 1 - «,(«) " . . 8 ■ 1 + Sr 2 ' W 



a x (s) 
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1 + 0/o 
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later sections 
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T 2 = T B I 





(44) 
(45) 

(46) 

(47) 

5.2 A Better Approximation jor the a Frequency Dependence in the Ebers- 
Moll Model 

Pritchard 41 has first suggested that a better approximation for the 
3-dB cut-off points of the a's or 0's are obtained if one inserts a factor 
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1.22 into the corresponding equations, i.e., 

a " = coshTAz W " " 1.22JC0 " ' (48) 

^out-of f measured 

This can readily be calculated from the fall-off behavior of the cosh 
expression while assuming jS^q 2$> 1. 

The same factor 1.22 appears in the corresponding expressions for 
<xj , (3 N and (3j . It is evident from (48) that this problem can be reduced 
to a matter of denning u atf . For less ideal transistors the factor is 
usually between 1 and 1.22. 

Higher-order approximations to the hyperbolic function commonly 
use two pole expressions or delay-producing excess phase terms. 

VI. THE CHARGE-CONTROL TRANSISTOR MODEL 

In analogy to the diode charge-control model we can establish a 
charge-control equivalent circuit for the transistor. To that end, we 
want to express all parameters in terms of the charge in the base. 

Three approaches appear feasible: A lumped Linvill model can be 
labeled in such a way that all elements appear as functions of charges 
rather than integrated charge densities of the form pAx . The two are 
proportional; the proportionality factors are of the form "electron 
charge times area". Most of the special circuit components of the 
Linvill model become current or voltage sources in the charge-control 
version. This procedure of converting a Linvill model into a charge 
control model can readily be applied to higher-order Linvill models. 

A second approach is to use the Ebers-Moll principle of superposition 
whereby two charge-control diode models plus the corresponding col- 
lecting currents can be joined to form the transistor model. This ap- 
proach is essentially limited to the first order of approximation. Two 
seemingly different, but fully equivalent and easily convertible models 
result. 

The third and classic approach to charge-control theory, originated 
by Beaufoy and Sparkes, 3 is basically mathematical. Through integra- 
tion of the continuity equation the carrier density as a variable is 
replaced by the total charge in the base. Certain simplifying assump- 
tions have to be made to obtain a relation between currents and charges. 
In essence, these assumptions are equivalent to the approximations 
implied in the first-order Linvill and Ebers-Moll models as well as in 
the first-order charge-control equivalent circuits to be described below. 
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Some equivalent circuits have been presented in the literature, but 
they were less rigorous than the circuits described below in the sense 
that they cannot be used as complete networks. Additional knowledge 
of the physics of the device is required to use these models. Extension 
to higher-order models in the Beaufoy-Sparkes approach is accom- 
plished through increased physical and mathematical complexity and 
not through more complex network topology as in the Linvill model 
or the charge-control model to be described. 

6.1 The ir-Version (Base-Controlled Version) of the Charge-Control 
Equivalent Circuit 

In the lumped Linvill 7r-model of Fig. 16, the base charge distribution 
is approximated by two levels of carrier density. This is illustrated 
in Fig. 18. p r (t) is constant over the length Ax = W x , and p 2 (t) is 
constant over the length W 2 , where W x + W a = basewidth W. The 
total charge in the two sections follows with (34a) and (34b) as 

qi (t) = Pl (f)WMi = PnoW 1 eA 1 [exp \\v.(t)} - 1] (49) 

<? 2 (/) = p 2 (t)W 2 eA 2 = p n0 W 2 eA 2 [exp {\v e (t)} ~ !]■ (50) 

Using the definitions of the elements given in Fig. 16, one can calculate 
from the Linvill model in Fig. 16 the currents through H el , H e2 , 
and H d and obtains 



*,(/) = H el pS) = PiW = ~T~ = ~ 

T\ T\ T B N 

i 5 (t) = H c2 p 2 (l) = p 2 {t) = —— = *- 

T 2 T 2 T B i 



(51) 
(52) 



*- p, = p(o) 




fc|*-P2=P(W) 



Fig. 18 — Excess carrier distribution in the transistor base as used or implied in all 
first-order transistor models, (a) General case pi ^ p(0), p 2 ?* p(W). (b) Commonly- 
used choice pi = p(0), p 2 = p(W). 
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!,(/) - H d [ Vl {t) - p 2 (t)] 

eAD 



W 



w»-*»-¥(*ft-ift)- w 



Using the more familiar Ebers-Moll notations and the relations found 
earlier in (42) through (47), i a can be expressed as 

m = — <is(t)-^ qi (t). (53b) 

Thus, the three current sources in the charge-control model are found 
and related to the Linvill model by means of (51) through (53). 

The remaining two branch currents i 2 and i A are obtained from 
Fig. 16 as 

m - U.W, ^f 1 = '^ (34) 

ii(t) = eAiWi *nM = <m. (55) 

These equations describe two stores S N and <S' 7 , whose properties have 
been described in Section 2.4. 

The conversion between the two models will be summarized and 
further discussed in Section G.4. 

The voltage sources for the junctions follow from (34), (51), and 
(52) as 

W «-In(l+^)-ln(l+-gM (56) 

MO = ln(l+^Uln(l+-^U. (57) 

With the help of (41) through (47) that link the constants used in 
the Ebers-Moll model to those in the Linvill Model, (56) and (57) 
can be rewritten as 

L Tjtti ' / E0 /(l — a,vott/n)J 

ImM = In ("l + 3M X - 7 / . 1 . + fa . 1 • (59) 

L T, n I co/ {I — Oi N na I0 )J 

(The reader may prefer to derive the constants directly from the 
steady-state Ebers-Moll model in (40) by considering the limiting 
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cases v t = and v c = 0.) With the addition of (58) and (59) the equiv- 
alent circuit in Fig. 19(a) is completely defined. 

It is customary and useful in charge-control work to define additional 
parameters t bn , t C n , t E i , t C i . We define their relationship as follows: 

(60) 
(61) 



(62) 





Tbn Ten 

o = = t C n 

HNO a N0 




Tbi t ci 

— = = TEI 

Pio oc I0 


Since, as usual, 






&N0 




am - 1 + to 


and 






010 
a '° - 1 + /3,n ' 


it also follows that 






Ten Tbn Tcn 


and 






_L = X + J_. 




Tci t bi t ei 



(63) 



(64) 



(65) 

(The classic definition of these time constants after Beaufoy-Sparkes 
will be discussed in Section 6.5.) The subscripts B, E, and C stand 
for base, emitter, and collector, respectively. The subscripts N and / 
on the time constants and on the charges have been chosen to indicate 
the normal and inverse transistor operation. Many authors 0,1 use 
F (forward) and R (reverse) instead of N and /. Since F and R are 
commonly reserved for diode forward and reverse currents, and since 
such currents can flow in each of the two junctions, the different nota- 
tions N and /, as proposed by Ebers and Moll, appear more appropriate. 

In Appendix B, the notations for the stored charges and the time 
constants used in this paper are related to those used in a recent book 
published by the Semiconductor Electronics Education Committee; 11 
they are also compared with the notations and definitions used by 
Beaufoy and Sparkes. 

The additional time constants do not add any additional degree 
of freedom. But it is advantageous to use "base" notations when 
controlling base current, i.e., in common-emitter or common-collector 
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Fig. 19 — Charge-control equivalent circuit for transistor in first-order approxi- 
mation, shown in two equivalent and convertible forms: (a) x-version, Linvill type, 
(b) T-version, Ebers-Moll type. 
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connection, and to use "emitter" and "collector" notations when 
emitter and collector forward currents are injected, such as in common- 
base connection. 

The model obtained in Fig. 19(a) maintains the most valuable 
property found in the Linvill model, namely the close relationship 
between the physical processes and the circuit elements. For example, 
T By and t bi are the recombination times on the emitter and collector 
side, respectively, and t cn and t ei are the diffusion time constants 
for the charges injected from the two junctions. Junctions and base 
are represented by individual sections within the equivalent circuit. 
This separation makes it easy to expand the model and to take other 
effects into account. 

6.2 The T-Version {Emitter-Collector Controlled Version) of the Charge- 
Control Equivalent Circuit 

In complete analogy with the derivation of the Ebers-Moll model 
in Fig. 17(a) we can take two diode charge-control models back to 
back and add current sources on the collector and emitter side, which 
are a N0 and a I0 times the diode currents. 

For the simulation of the junctions, we are left with two alternatives: 
One is to convert the corresponding expressions in the Ebers-Moll 
model in Fig. 17(a) into charge functions; the other is to use the expres- 
sions in the charge-control 7r-model (which are equivalent to the Linvill 
model), but replace the /3-notations by a-notations according to (GO) 
through (63). The first-mentioned alternative for simulating the voltage 
sources would amount to simply substituting the diodes from Fig. 17(a) 
for the voltage sources in Fig. 19(b). The property of charge control 
would not be simulated. The second procedure is therefore chosen; 
it yields 

■•«> = k ln (* + ~ j /n {t) —) W 

X \ tea-iw(1 — at,aan)/ 
v c (t) = \ In (l + T „? (0 r) ■ (67) 

X \ T C jico/(l — a.vott/n)/ 

Thus, the equivalent circuit in Fig. 19(b) is obtained. 

As far as the current relations in the models are concerned, the 
main difference between the charge-control T-model and the Ebers-Moll 
model is that the frequency dependence is simulated by a mathematical 
expression in the Ebers-Moll model, and by an additional network 
branch in the charge-control model. This is analogous to the option 
existing in small signal models where one can represent the frequency 
dependence either with an appropriate RC circuit, holding a frequency 
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independent, or alternatively, with a frequency dependent a in the 
collecting current source. 

The equivalency of the charge-control model with the Ebers-Moll 
model exists only for the relations between the currents. It can be 
shown readily that the following relations must be satisfied to establish 
equivalency: 
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(69) 
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Till - — 

0)gi 03 a I 


(71) 



All o/s must be replaced in these equations by the corresponding 
«/ 1.22 if the w's correspond to the measured 3-dB gain fall-off points, 
and if the better approximation mentioned in Section 5.2 is to be 
included in the Ebers-Moll model, provided the particular transistor 
follows the underlying theory well enough. 

G.3 Conversion Between the Two Proposed Charge-Control Models 

The identity between the two charge-control models, presented in 
Figs. 19(a) and (b) can best be proven by converting one model into 
the other. 

To convert the x-model into the T'-model one first adds a branch 
current i d both into and out of the base point B' and splits i d up into 
its two components. The resulting circuit diagram is shown in Fig. 20. 
The two parallel current sources proportional to q N (t) on the left side 
can then be combined into one current source. Likewise, the two current 
sources proportional to q,{t) on the right side can be combined. If 
with the help of (60) through (63), one now relabels all current sources 
in terms of a N and a, instead of fi N and /3, and extends the upper current 
sources beyond the voltage sources, one obtains the model in Fig. 19(b). 

6.4 Summary of the Conversion Equations between the Linvill and the 
Charge-Control Model 

6.4.1 Conversion Equations for the First-Order Transistor Model 

In (49) through (59), the charge-control r-model was derived from 
the Linvill model. With the help of the defining equations for the 
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Fig. 20 — Intermediate step used in the conversion from the ir to the T charge- 
control model, demonstrating equivalency between these two models. 

Linvill model elements, the constants in the charge-control model can 
be calculated as a function of the Linvill combinances, storances, and 
diff usances. For the relation between the Linvill 7r-model of Fig. 16 
and the charge-control 7r-model of Fig. 19(a), such calculations yield 



tbn = 



TBI = 



TCX 






§1 



1 E0 — PNO 



H cl H c - + (Hci + H cl )H d 



_ H fi H el + (//,. + H ea )H d 



(72) 
(73) 
(74) 
(75) 
(76) 
(77) 



Note that (72) through (75) reveal that the five parameters in 
the Linvill model lead to only four parameters in the charge-control 
model. The one degree of freedom that is lost in the charge-control 
model is the conversion factor from current to carrier density; con- 
version of the charge-control model into a Linvill model is only possible, 
if one of the five Linvill parameters is known. This is tantamount to 
saying that one needs some information on the geometry of the device 
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such as the value of one, or in low-gain units, both of the two base 
volume sections AiWi and A 2 W 2 • 

6.4.2 Conversion Between the Linvill and the Charge-Control Model for 
an Arbitrary Number of Base Sections 

In higher-order approximations for diodes or transistors, the param- 
eters of the Linvill and the charge-control models, as defined in Fig. 21, 
are related by the equations 
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Fig. 21 — A Linvill and a charge-control equivalent circuit for a junction and part 
of a multisectional n-region, with indication of the notations used in converting one 
model into the other. 
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6.5 The Transistor in Saturation 

In all lumped transistor models (Linvill, Ebers-Moll, or charge- 
control type) the charge in the base is explicitly or implicitly broken 
up into the charge q N injected from the emitter under normal operation 
and the charge q r injected from the collector under inverse operation, e.g., 
in saturation. This was illustrated in Fig. 18. 

When the transistor is overdriven into saturation with a base current 
larger than I c ..t/ftro , the two stores q„ and q r do not change by exactly 
equal amounts, i.e., 

dq N _ dqi 



where, by definition, 



'h> 



di, 



= %R — 



ft 



(83) 



This is illustrated in Fig. 22. It can be calculated from any of the two 
models of Fig. 19 that, under steady-state conditions, the excess charges 
in the two stores are related to the excess base current by the expressions 



Mh = Qr = 



1 — O-vnttM 



/«„,,.,, 



(84) 



AQ N = 



TKN 



/„™ H . 



(85) 




Fig. 22 — The transistor in saturation, (a) Actual distribution of excess minority 
carriers, (b) Lumped approximation. The tt and the T models use q N and qi with 
lifetimes t bn and t bi ; the Beaufoy-Sparkes model uses Q N sat and q BS = Ag A r + q, 
with lifetimes t bn and t s , where t s is as defined in (88). 
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From this it follows that 

AQ,v t en Ton a 



Qi otf/oTcr (XiqTei at m io aty 



^ 1. (86) 



Since a ri) < 1, the charge-up ratio is somewhat larger than the ratio 
of the diffusion times of the normal and inverse transistor. 

The rate at which the two stores charge and discharge in saturation 
because of external step disturbances is described by the eigenfunction 
of the system 



s 2 + s 



1 + § m 1 + ft 



TBN T bi _ T BN T H i 



+ 1 + &a + gffl . o. (87) 



If &nqTm/t bn >>> 1, the two poles arc far apart in frequency. Further- 
more, the high frequency pole contributes in most nonoscillatory cases 
little to the overall response. The higher pole or, alternatively, the s 2 
term can then be neglected and a single time constant results described 
by 

_ (1 + /3.vo)tb/ + (1 + fto)T B .y Tg.y + Tci (QO s 

1 + P.vo + Pio 1 — oc. vo a IO 

Using w-notation, one obtains the expression given by Ebers and Moll 

rs = "■; + "-' • (88b) 

For largo p\ v and small /3/ , r. s - is approximately equal to 

ts^tJi + — )• (88c) 

If t E h « 7-f/ , i.e., if the carriers diffuse more easily from the emitter 
to the collector than vice versa, then the recombination rate t ui on 
the collector side is mainly responsible for the overall decay of the 
excess base charge. 

6.5.1 Storage Time Calculations 

For first-order storage time calculations with the transistor driven 
into a steady-state saturation condition by means of an excess base 
current I „ exceB8 , one can simplify the charge-control model to the one 
shown in Fig. 23. Storage time is the time it takes to deplete the store 
which is charged to a value of 

QBS - Qr + AQ* = Z,«.„.r, = I,„... ^ + "^ (89) 

1 — a,v„a/o 
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Fig. 23 — Single-pole equivalent circuit for saturated transistor after Beaufoy- 
Sparkes. 

while, at the same time, this charge is exposed to an effective recombina- 
tion time of t s as given in (88). In this form the model is identical 
with the classic Beaufoy-Sparkes model for the saturated transistor. 
In the general case one must refer to the complete model. 

6.6 The Beaufoy-Sparkes Charge-Control Model 

In the classic approach to charge-control theory, the starting point 
is, like in the diode case, the integration of the continuity equation (8). 
In comparison with the integration performed for the diode in Section 
2.3.1, the upper limit of integration has to be changed to x = W. The ex- 
pression 

*,«>, - i v (W, t) = -^- + — 

obtained from the integration becomes that for the base current under 
normal, nonsaturated operation: 



i .f. = dq N (t) q N {t) 
at Tbn 



(90) 



q N (t) is the total charge in the base. The next step being made is again 
the approximative assumption that the carriers redistribute themselves 
so quickly, that we can always assume steady-state distribution. 
(See also Section 2.4.) Mathematically, this means that in normal 
transistor operation both p(0,t) and i (t) are proportional to the base 
charge q N (t). It can be seen from Fig. 19 that the same assumption 
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is implied in the two charge-control models presented there, despite 
the fact that they were derived through entirely different procedures. 
(Instantaneous redistribution is, however, not implied in models that 
use more than one ir or T section for representing the base.) 

The time constants are denned in the classic charge-control theory 
on the basis of the above-mentioned assumption of instantaneous 
carrier redistribution, i.e., in steady state 

Q.v Q.v Qy /n-i„\ 

t B n = J — , t C n = 7 — , r EN = -. — (91a) 

■* B.V * CN * EX 

and dynamically 

Q.v , "<7.v • Q.v ■ • i • /mu 

lB!f = - J — + -JT , ICN = , tEN = iRN + ICN • (91b) 

TbN «*■ Tcn 

The remaining three time constants can be defined likewise for the 
inverse transistor. Narud, et al 9 have used such definitions in an equiv- 
alent circuit for the charge-to-current relations in the transistor. Beaufoy 
and Sparkes discussed this possibility in their original paper 3 but chose 
to present two separate charge-control models, one for the normal 
active operation and one for saturation. In normal operation, the 
charge q N called "q B " is bounded by the value reached at the edge 
of saturation : 

Qb ^ Qy ,at , where Q, v sat = -J 854 r BN . 

Pm 

In their saturated model, all excess charge which exceeds Q N sxr is 
lumped into one store rather than two; this charge "q n s" has a lifetime 
ts = 9bsAb exec, which is identical with t s as defined in (88). 

By lumping Aq N = q N — Q N SAT and q r into q B s , the Beaufoy and 
Sparkes arrangement provides only a minor short cut for calculating 
storage time, while sacrificing not only some of the physical under- 
standing, but also the possibility of mutual conversion with the other 
models. \ii relations have been given that would express the junction 
voltages in terms of the charges in the stores, and recourse must be 
taken to the Boltzmann equation to find expressions for the voltages. 

Throughout the literature the charge-control concept has been used 
primarily as a mathematical-physical tool. Extensions to higher-order 
effects are usually made by improving the simple continuity and 
transport equations stated in (8) and (9) and then carrying out the 
corresponding integration for the specific application. 
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VII. SOME REMARKS ABOUT THE EQUIVALENT CIRCUIT TYPE CHARGE- 
CONTROL APPROACH 

7.1 Use of the Charge-Control Models 

It is believed that the first-order approximation to a charge-control 
model in the form presented for the transistor in Fig. 19, combines 
the main advantages of the three basic approaches to modeling. The x 
and the T-models are as easy to handle from an equivalent circuit 
point of view as the Ebers-Moll model. Instead of frequency dependent 
a's and /3's, one additional current branch exists for each side of the 
transistor. Circuit problems are solved in the usual way by means 
of loop and node equations. The charges q s and q { appear as circuit 
parameters which can either be calculated, if so desired, or else, elim- 
inated in the algebraic process. The store elements in the circuit are 
clearly defined by the circuit properties given in (22). 

The model provides all the features that have made the charge- 
control concept attractive in the past: quick estimates of switching 
times by integrating the base current and equating with the charges 
needed to fill and deplete the stores. The general base current equations 
of charge control are directly read from Fig. 19(a) as 

q» , dq N , q, , dq, C T s &U _l C ™ dv < . (QO\ 

lB = 7^ + ~dT + 77, + ~dJ + dt + dt m 

Of course, there is no restriction to step inputs. The chore of cal- 
culating responses to a nonstop input is transformed through the model 
into a circuit problem. In complex cases the help of a computer will 
be required. 

Due to its direct relationship to the Linvill model, the charge-control 
model lends itself quite readily to extensions based on the physics 
of the device. This will be discussed in Section 7.3. 

7.2 Piecewise Linear Approximation of the Logarithmic Voltage Function 
The logarithmic voltage functions for the junctions are of the form 

i + Tur^ — *T (93) 

For most practical cases, this can be approximated by piecewise linear 
functions, like in the diode case of Fig. 10. Except for small values 
of q/r, i.e., q/r not »/ /(l - a m a TQ ) , one obtains 

dv = -I. (94) 

dlq/r) \q 
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Thus, the slope can be represented by a resistor r/\q, which may, 
like in Section 2.4.4, be taken as the average value 

1 (95) 



where I max is the maximum forward junction current. 

It can be shown that in models which use the exponential relation- 
ship, the expressions of the form ln(l + x) can be replaced by just 
In z if one simulates the majority carrier currents by special current 
sources as follows: 



and 



1 ~ a in 
1 — a.voa/o 

1 — tt.vo 
1 — a.voa/o 



I co (~/c<>) from internal base to collector 
I eo (~0) from internal base to emitter. 



This transformation is rigorous only at dc. However, in a piecewise 
linear analysis, as discussed above, the addition of one current source, 
namely I co , becomes mandatory if the model is to be valid at very 
small collector currents. 

7.3 Extensions of the Model 

7.3.1 Path Impedances, Leakage Resistors 

Like in the Linvill model, junctions and base material are clearly 
separated in the charge-control model. Therefore, it is a straight- 
forward procedure to add series path resistors, series inductances, or 
leakage resistances to models like the ones in Figs. 9 or 19. 

7.3.2 Junction Capacitors 

It has been indicated by the dashed lines in Figs. 9 and 19 how the 
junction capacitances are to be incorporated into the model. They are 
properties of the junction, but their currents flow as majority carrier 
currents through the bulk material. Hence, in Fig. 9, for example, they 
must be connected across the whole n-region. (Connecting directly 
across the voltage source would have no effect on the external prop- 
erties.) In Fig. 19 they lead to the internal base point. 

7.3.3 Higher Order than ir-Transislor Models 

Another desirable expansion may be to replace the ir structure of 
Fig. 19(a) by a double ir or by some other higher-order approximation 
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to the original differential "transmission line". This is of special im- 
portance, if emphasis is to be placed on charge redistributions in the 
base. By extending the model in such a way, the restricting assumption 
of instantaneous carrier redistribution is no longer implied. A qualitative 
example of an elaborate planar or mesa transistor model is given in 
Fig. 24. 
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Fig. 24 — Example of an elaborate high-frequency planar or mesa transistor 
equivalent circuit. 

7.3.4 Drift Fields 

If the charge-control model under consideration is being developed 
on the basis of physical phenomena such as in the model of Fig. 19(a), 
the contribution from drift effects may be represented in the same way 
as has been proposed by Linvill [Ref. 7, Sections 2, 3]. As a direct con- 
sequence of the transport equation (2), drift can be represented by a 
current source added in parallel to the diffusion current source. In terms 
of the r-g-c transmission line representation, discussed earlier, drift con- 
sideration amounts to a resistor in parallel with the series diffusance 
resistor r. This was used in a recent paper by Bloodworth. 26 

Alternative methods of representing drift effects in conjunction with 
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conductivity modulation are presently being investigated; results will 
be published later. 



7.3.5 Base Width Modulation 

Base width modulation can be taken into account by replacing the 
basewidth, especially the collector section W 2 , by an expression 
W 2 (l + A), where A is some function of the junction voltage. Equa- 
tions (52) and (53) show the dependence of the branch currents on W 2 , 
from which we can readily derive the required modification of the 
charge-control model in Fig. 19(a). 

7.3.6 Multiple-Layer Devices, Multiple Storage 

In accordance with Linvill's proposal, storage in more than one 
region can be simulated by considering that the minority carrier cur- 
rent on one side of a junction becomes the majority carrier current 
in the adjacent region. An example is shown in Fig. 25. This figure 
represents the charge-control model for an npnp device. Avalanche 
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Fig. 25 — Charge-control model for npnp device. The model for an npn-transistor 
with storage in the collector can be obtained from this model by omitting the part 
to the right of the dashed or the dotted line. 
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multiplication may be considered by adding a multiplication factor M 
to all hole and electron currents flowing through the junction of interest 
(usually the center junction J 2 ), as indicated in Fig. 25. M is a function 
of the voltage across the junction. 

By omitting the last electrode, an npn transistor with charge storage 
in the collector is obtained. 

7.4 Establishment of a Large-Signal Model 

The question naturally arises as to how one arrives at a numerical 
model. There is no clear-cut answer to this question, since the procedure 
to be taken depends on whether the informations available are pre- 
dominantly physical or electrical in nature, whether a computer is 
available or not, etc. The following outline can, therefore, only be 
considered as a typical example. 

(i) Obtain dc measurements which yield information on junction 
characteristics and electrode resistances. All measurements must be 
made under widely differing drive and load conditions. 

(ii) Add information from device manufacturer to establish first- 
order dc model. (If necessary, convert to Linvill model.) 

(Hi) Add dynamic parameters, such as capacitances, as far as they 
are known and establish first-order dynamic model. 

(iv) Use computer to improve numerical parameter values by 
matching frequency response curves or switching data in the active 
region with the model. 

(v) Use computer to match large-signal nonlinear switching data. 

(vi) Check model with switching measurements under different con- 
ditions, such as extremely low, extremely high and medium input and 
output impedance levels for various drive conditions. Improve model 
basically and numerically as necessary. 

For purposes of device design, more emphasis is generally placed on 
the simulation of higher-order effects than in model building for circuit 
design where, especially in the case of integrated circuits, it is necessary 
to trade accuracy for simplicity. 

VIII. CONCLUSIONS 

The differential Linvill model stands out among all models as the 
most perfect one. Whereas the lumped Linvill model is the most suitable 
model for the device physicist, the circuit engineer usually prefers a 
more circuit oriented approach. It is felt that the charge-control equiv- 
alent circuit approach is well suited to combine the main advantages 
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of the various models: It is as easy to handle as the Ebers-Moll model, 
yet bears the close relationship to the physical phenomena of the 
device inherent in the Linvill model. It can also be extended easily 
to include higher-order physical phenomena. 

Despite the difference in basic philosophy underlying the creation 
of each of the three classic modeling concepts (such as lumping, super- 
position, and integration), they are equivalent with respect to their 
current relations and to all dc properties. When compared at the same 
level of complexity, equivalency with respect to time dependency of 
the junction voltages exists between the two charge-control models 
and the Linvill model, but not between these models and the Ebers- 
Moll model. 

In the Ebers-Moll model the effect which storage exercises on voltage 
cannot be included. The hydrib use of both time and frequency domains 
in the model may also be felt as a disadvantage in some applications. 

At the first-order level of approximation, the charge-control equiv- 
alent circuit can be converted into the Ebers-Moll model, the Beaufoy- 
Sparkes model, and into the Linvill model (in the latter case the base 
volume is a constant which must also be known). Thus, the charge- 
control model serves as a bridge between the various models. This can 
be very useful in establishing a model, since both physical and electrical 
information can be incorporated easily into the model. 

The diode charge-control model has been found very useful for 
analyzing storage effects in diodes. 

Because of the close relationship to the physical phenomena in the 
device, extensions to larger complexity can readily be accomplished. 
We may interpret the charge control equivalent circuit as simply a 
circuit-oriented form of the Linvill model. The basic ideas and pro- 
cedures that are used in converting diode and transistor linear models 
into equivalent charge-control models can be applied to many other 
semiconductor devices. 

APPENDIX A 

Switching Time Calculations for Ideal Charge-Storage-Step-Recovery 
Diodes 

(Example for use of charge-control model) 

A.l Equivalent Circuit (See Fig. 26) 

A. 2 Generator Source Current (See Fig. 27) 

A.3 Diode Model (See Fig. 28) 
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Fig. 26 — Equivalent circuit for charge-control model. 
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Fig. 27 — Generator source current. 
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Fig. 28 — Diode model. 



A. 4 Forward Operation 



From this follows 



m = m + sQ(s) = u. 

T S 



Q(s) = 



ia + i 



q(t) = rl,[l ~ exp(-//r)] 
Q(Q = rI F [l ~ exp(-/ P /r)] 
A. 5 Reverse Operation 

For simplicity of writing, t = t p will now be referred to as t = 0: 
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1(8) - m + sQ(s) - Q(Q = - ! -f- 



From this follows 



Q® = 



SQ(Q - In 



is + \ 



q(0 = Q(Q exp i-t/r) - tI b [1 - exp (-t/r)]. 
Step recovery occurs at t = T x , when q = 

exp (-T x /r)[Q(Q + r/„] = tI r 

T z = r In TlH- ^1 = r In (l + 7^ [1 - exp (-/ p /r)] 

A.6 Graphical Representation (See Fig. 29) 

q(t) 

T 
I . 

Q(t p ) 




-Irt 



Fig. 29 — Graphical representation. 



A. 7 Special Cases 

(i) t, -» 00 : Q(0 - i> 

r.-ri»(i + g 



(it) I R » / F : r, = 



QUA 

Ik 



(Hi) I R » If and t„ -> » : T m = -f- r. 



(90) 



(97) 
(98) 
(99) 
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APPENDIX B 

Comparison o] Notations 

Table I lists comparisons of the notations used in this article. The 
first column lists the notations used in this article while column A 
lists those used by Beaufoy and Sparkes. 3 Column B lists the notations 
used in Physical Electronics and Circuit Models by P. E. Gray, et al; 11 
SEEC Series, 2. 

Table I — Comparison of Notations 



This paper 


A 

Beaufoy-Sparkea 


B 

SEEC 


ft* 


q„(for q x 


< Qnsat 


only) 




'//•• 


<ll 
j Aq N + q r ~\ 












'nt 


1 where Aq N = q lW - Q NS at) 






qBs 








tbn 






T B 






Till- 


T CN 






T c 






Tp 


TEN 






T B 




7 


\TBF Tp] 


TBI 












TBR 


TCI 










7 


7-L + f) 

\TBR Tr) 


TBI 












TR 


T.S 






T s 






TSL 



LIST OF SYMBOLS 



Lower-case letters are used for time variables, capital letters are used 
for steady-state values or Laplace transforms of values. 



A, A, , A 2 

■^•11 > -^12 > -^21 

a, b, c, K 
c; c 
tfu C/j»« 

D p 

D n 

e 

E 

Q',9 



cross-sectional areas 

four-pole parameters 

constants 

analog capacitance; same per unit length 

emitter and collector junction capacitance, 

respectively 

hole diffusion constant 

electron diffusion constant 

magnitude of electronic charge 

electric field intensity 

analog shunt conductance; same per unit 

length 
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H_ c , H el , H c2 lumped combinances 

H c combinance per unit length 

II, i lumped diffusance 

1/Hj reciprocal of diffusance, per unit length 

h , i>i ) ia , U , io network branch currents 

i B base current 

i>a excess , Ib exees., excess base current in saturation 

i c , I c collector current 

ig , I B emitter current 

I P , I R forward and reverse diode switching current, 

respectively 

i„ electron current 

i v hole current 

I s diode saturation current 

ic , Ic y is j Ie 1 network branch currents as defined in 

ice , Icp , i-EF , Ief /Figs. 17(a) and 17(b) 

I , I co , I Eo , Icfo , Iefq dc junction saturation currents 
I c Sat collector current in saturation 

i Hc , I IU , i s currents through combinance, diffusance 

and storance, respectively 

Ibn , isN ) icN base, emitter, and collector current in normal 

transistor operation 

j n electron current density 

j p hole current density 

Ic Boltzmann constant 

L„ diffusion length for holes 
m, nii , m 2 constants relating lumped carrier density 

to carrier density at junction boundary 

n electron density; excess electron density 

n„ excess electron density in p-region 

n v , n p0 values of n and n p in thermal equilibrium 

p, P(s) hole density or excess hole density 

p n excess hole density in n-region 

p , p„ value of p and p n in thermal equilibrium 

q, Q charge 

Qi } Q2 , Q m lumped charges in base region 

q N , Q N charge in normal store 

q r , Q T charge in inverse store 
Ag\ v , &Qn , Ag, , AQi , Qbs additional charges stored due to saturation 

Q, v BAT limiting value of Q. v , reached at edge of 

saturation 
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total minority carrier charge in equilibrium 
analog resistor, same per unit length 
small signal junction resistance 
Laplace operator 
stores = storances 
storance per unit length 
time 

absolute temperature 
voltage 

externally applied junction voltage (ex- 
cluding resistive drops) 
collector and emitter junction voltages 
lengths denoting sections in neutral region 
base width 

neutral region length variable 
characteristic impedance 
normal and inverse ac current gain in 
common-base connection 
dc values of a N and a, 
normal and inverse ac current gain in 
common-emitter connection 
dc values of pV and /3 7 
transmission line propagation constant 
short for e/kT 

electron and hole mobility, respectively 
recombination time constant in p-region 
diffusion time constants 
approximative effective recombination time 
constant for excess charge in saturation 
recombination time in base under normal 
operation 

recombination time in base under inverse 
operation 

normal collector and inverse emitter diffu- 
sion time constants, respectively 
normal emitter time constant (= l/co aN ) and 
inverse collector time constant (= 1 /«„/), 
respectively 

common-base angular cut-off frequencies 
common-emitter angular cut-off frequencies 



Qio, 


q m o 


r; f 




r t 




s 




S,S 

s 


i , S 2 , S m 


t, t x 


,T 


T 




v,V 




»«t 




v e , v„ 


w 1 


w 2 


W 




X 




Z 




oc N , 


CLl 


OCKO 


«70 


frv, 


ft 


PiVO 


fro 


7 




X 




fin, 


Mp 


T, T„ 




K, 


T* , T b , T12 


T, 




T* = 


= Tbn 


T 2 = 


= TBI 


Ton 


, Tei 


ten 


, T C J 


W„iV 


, w„, 


W/JA, 


o>m 
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